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In this article we study in detail the prospects of determining the infrared finite QCD effective 
charge from a special kinematic limit of the vertex function corresponding to three background 



gluons. This particular Green's function satisfies a QED-like Ward identity, relating it to the gluon 
D 

propagator, with no reference to the ghost sector. Consequently, its longitudinal form factors 

may be expressed entirely in terms of the corresponding gluon wave function, whose inverse is 

> 

"^j- ■ proportional to the effective charge. After reviewing certain important theoretical properties, 

ON. 

1/"} , we consider a typical lattice quantity involving this vertex, and derive its exact dependence on 

(N 

T-f- ' the various form factors, for arbitrary momenta. We then focus on the particular momentum 

O : 

configuration that eliminates any dependence on the (unknown) transverse form factors, projecting 

out only the desired quantity. A preliminary numerical analysis indicates that the effective charge 

S^ is relatively insensitive to the numerical uncertainties that may afflict future simulations of the 

h' 

aforementioned lattice quantity. The numerical difficulties associated with a parallel determination 
of the dynamical gluon mass are briefly discussed. 

PACS numbers: 12.38.Aw, 12.38.Lg, 14.70.Dj 



I. INTRODUCTION 



The ongoing quest for a deeper understanding of the nonperturbative features of the basic 
Green's functions of QCD has been benefited considerably by the constructive interaction 
between lattice simulations and the Schwinger- Dyson equations (SDEs). One of the central 
issues in this search is the systematic study of how the effects of dynamical mass generation 
manifest themselves in the nonperturbative structure of the propagators of the fundamental 
fields (quarks, gluons, and ghosts). In particular, the concept of a momentum-dependent 
gluon mass |l| nas received renewed attention, because it provides a natural explanation 
for the infrared finiteness of the gluon propagator, A(g 2 ), and ghost dressing function, 
F(q 2 ), observed in large- volume lattice simulations carried out in the Landau gauge, both 
in SU(2) [2] and in SU(3) [3|. 

Especially interesting in this context is the definition and properties of the QCD effective 
charge, as well as its interplay with the gluon mass. In fact, as has been explained in detail 
in a series of works |4j-|6|, if the finiteness of the aforementioned Green's functions is to 
be attributed to the dynamical generation of a gluon mass [i.e., A _1 (0) = m 2 (0)], then the 
corresponding effective charge, when properly defined, saturates also at a finite, nonvanishing 
value |l|, |5j. Within this picture, the so-called "freezing" of the charge is a consequence of 
the presence of the gluon mass in the corresponding logarithms, which are cured from the 
perturbative Landau pole, are well-defined for all physical momenta, and acquire a finite 
value at q 2 = 0. It is important to emphasize that this characteristic property of the QCD 
coupling, which endows the theory with an infrared fixed point, has been advocated by a 
large number of very distinct theoretical and phenomenological approaches, see |ll. la. I7H14| 
and references therein. 

A self-consistent framework for the study of the effective charge is provided by the fusion 
of the pinch technique (PT) |ll. Il5l-ll9| with the background field method (BFM) |20(, known 
in the literature as the "PT-BFM" scheme 2jJ-|23|]. The natural starting point in this context 
is the PT gluon propagator, A(g 2 ), which is known to coincide with the two-point function 
of two background gluons. Note that A(g 2 ) and A(g 2 ), as well as their individual kinematic 
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25| [see, e.g., Eq. (1279]) 



components introduced below, are related by powerful identities 

involving an auxiliary two-point function, which, in the deep infrared, practically coincides 

with the ghost dressing function [26J , |27J • 



In the presence of a dynamically generated mass, the (Euclidean) A(q 2 ) assumes the 
form A _1 (g 2 ) = q 2 J(q 2 ) + rh 2 (q 2 ), where the first term corresponds to the "kinetic term", or 
"wave function" contribution, while the second denotes the momentum-dependent mass [28 ]. 
Then, the (infrared finite) effective charge a(q 2 ) may be defined as a(q 2 ) = a s (n 2 )J~ 1 (q 2 , /i 2 ), 
where /i is the renormalization (subtraction) point chosen, within an appropriate renor- 
malization scheme [28j. By virtue of the QED-like Ward identities (Wis) satisfied in the 
PT-BFM framework, this latter quantity is formally renormalization group (RG)-invariant 
(//-independent). Note that the freezing property of this a(q 2 ) hinges crucially on the fact 
that the infrared finiteness of A(g 2 ) has been accounted for by the inclusion of the fh 2 (q 2 ) 
in the above decomposition; instead, the alternative definition in terms of the gluon dress- 
ing function Z(q 2 ) = q 2 A(q 2 ), namely a(q 2 ) = a s (fi 2 )Z(q 2 , fi 2 ), gives rise to a RG-invariant 
effective charge that vanishes trivially at the origin [5]. 

Evidently, in order to proceed further with the determination of the physical effective 
charge, one needs information on the behavior of J(q 2 )- To be sure, J(q 2 ) could be obtained 
from Eq. (J2.13P if its conventional counterpart, J{q 2 ), were known; however, unlike F(q 2 ), 
this latter quantity may not be directly extracted from existing lattice simulations. Indeed, 
the lattice determines the full propagator A(g 2 ), but offers no direct information on the 
momentum dependence of the individual components. 

Similarly, within the nonperturbative PT-BFM framework, both J(q 2 ) and m 2 (q 2 ) satisfy 
independent (but coupled) non-linear integral equations, which are obtained from the SDE 
of the A(g 2 ) following a well-defined procedure, developed in a series of recent works [29|, 



30j. These two equations may, in principle, determine the complete dynamics of these two 
quantities. In practice, however, one is considerably limited by the fact that the main 
ingredients entering in them are the (largely unexplored) fully dressed three- and four-gluon 
vertices, for arbitrary values of their momenta, and, as a result, only approximate solutions 
may be obtained. 

It would be clearly desirable, therefore, to determine J{q 2 ) from an approach that is 
ab-initio exact, in the sense that it does not involve any field-theoretic approximations. To 
that end, in this article we explore the possibility of extracting this important quantity 



from a possible lattice simulation of the PT-BFM three-gluon vertex 15ll3l|. to be denoted 
by T, defined as the one-particle irreducible part of the correlation function involving three 
background gluons (the prospects for a lattice nonperturbative formulation of the BFM have 



been recently revitalized due to the results of [32J, |33|). 

The fundamental reason why this vertex can furnish clean information on J(q 2 ) is the 



Abelian WI that it satisfies [15J. Specifically, to all orders in perturbation theory this WI 
involves only the difference of A _1 (g 2 ) = q 2 J{q 2 ), defined at the appropriate momenta; this 
is to be contrasted with the usual Slavnov- Taylor identity (STI) satisfied by the conventional 
three-gluon vertex [34], which involves, in addition, contributions from the ghost-sector of 
the theory, and especially from the so-called ghost-gluon kernel [351 ] . 

Of course, to properly account for mass generation, T must be supplemented by a spe- 
cial nonperturbative vertex, denoted by V, which contains the necessary poles to enforce 
gauge invariance in the presence of a gluon mass [29]. This particular vertex is completely 
longitudinally coupled, and its divergence furnishes precisely the missing mass terms that 
convert q 2 J{q 2 ) into a massive A _1 (g 2 ), thus maintaining the form of the original WI intact. 
However, when the full vertex r + V is contracted by three polarization tensors, as happens 



typically in lattice calculations in the Landau gauge [36], any reference on the (completely 



longitudinal) V disappears 29], and only the dependence on the nonperturbative J, defined 
at different momenta scales, survives. Then, a special kinematic limit, frequently employed 
in the lattice studies of three-point functions [36], converts the lattice quantity of interest 
into a function of a single variable, given simply by the first derivative of q 2 J{q 2 )- Finally, 
a straightforward integration of the lattice result over the relevant momentum interval fur- 
nishes J{q 2 ), and consequently a(g 2 ); this constitutes the central result of the present work. 

It is clear that the usefulness of the aforementioned exact result must be assessed within 
the context of a realistic lattice simulation, taking into account, to some extent, the practical 
limitations associated with such an endeavor. In particular, it is important to provide a rough 
estimate of the errors that the various numerical uncertainties may introduce to the effective 
charge and, subsequently, the gluon mass. A simple modeling of these effects reveals that 
the predictions obtained for the effective charge are rather robust, and that the induced 
deviations do not alter significantly its theoretically expected behavior. Instead, with the 
exception of the deep infrared, the extraction of the gluon mass is afflicted by important 
qualitative discrepancies. 

The article is organized as follows. In Sec. [TT] we define the basic quantities appearing 
in this problem, and summarize some important relations, characteristic to the PT-BFM 
framework. In Sec. lIHI we introduce the RG- invariant definition of the QCD effective charge 



and of the gluon mass; in particular, we put forth an interesting analogy between this 
latter RG-invariant mass and the standard constituent quark mass. In Sec. IIVI we present 
the relevant three-gluon vertex, together with its basic properties. Particular attention is 
payed to the linear (ghost-free) WI satisfied by this vertex, which allows for the complete 
determination of its longitudinal form factors in terms of the gluon propagator, with no 
reference to ghost Green's functions. Sec. |V] contains the main results of the paper. First, 
we briefly review the general prospects of simulating PT-BFM Green's functions on the 
lattice. Then, assuming that such a simulation can be actually carried out for the vertex, 
we show how the effective charge and gluon mass can be directly reconstructed from the 
data obtained in a commonly used kinematical limit. In addition, we carry out a detailed 
numerical study on how the unavoidable numerical errors of a possible simulation propagate 
to the relevant physical quantities. Our analysis shows that, while a faithful approximation 
of the effective charge can be generally obtained, the reconstruction of the gluon running 
mass is much more subtle, displaying large fluctuations due to the inevitable distortion of 
delicate numerical cancellations. Finally, our conclusions are presented in Sec. I VI I 

II. GENERAL CONSIDERATIONS 

The full gluon propagator iA"^(q) = 5 ab Afj, v (q) in the Landau gauge is defined as 

\„(q) = -zP^(g)A(g 2 ), (2.1) 

where 

P,M) = 9»»- q -^f ( 2 - 2 ) 

is the usual transverse projector, and the scalar cofactor A(g 2 ) is related to the (all-order) 
gluon self-energy H^ v {q) = P^ v (q)U(q 2 ) through 

A-V)=g 2 + *II(g 2 ). (2.3) 

It is advantageous to introduce the inverse of the gluon dressing function, J(q 2 ), defined 
as [35| 

A" V) = g 2 J(g 2 ). (2.4) 

At tree- level, J{q 2 ) = 1. Perturbatively, at one-loop, it is given by 

2 13C A g 2 (q 2 \ 

J(g) = 1 + ^6-r ln Uj' (2 - 5) 



where Ca denotes the Casimir eigenvalue of the adjoint representation [Ca = N for SU(N)], 
and the renormalization has been carried out in the momentum-subtraction (MOM) scheme. 
Evidently, J~ x {q 2 ) displays a Landau pole at q 2 = p 2 exp(—96iT 2 / lSC^g 2 ) . 

The generation of a dynamical gluon mass leads to the infrared fmiteness of the gluon 
propagator, to be denoted by A m (q 2 ). In particular, in Minkowski space one has 

A^(q 2 ) = q 2 J m (q 2 )-m 2 (q 2 ), (2.6) 

with m 2 (0) 7^ 0. The subscript "m" in J m indicates that the resulting expressions are 
regulated by the presence of m 2 (q 2 ). Specifically, after gluon mass generation, the J m (q 2 ) 
may be qualitatively described by 

2 13C A g 2 ( q * + p m \ q *) \ 
J.M ) = 1 + -^-r In [ ~ 2 ) , (2.7) 

with g 2 , the constant p, and m 2 (q 2 ) such that J m (q 2 ) > for all values of q 2 , thus avoiding 
completely the appearance of a Landau pole. 

The detailed dynamics that govern J m {q 2 ) and m 2 (q 2 ) are determined by two coupled 
integral equations, obtained from the SDE of the gluon propagator. Specifically, after a 
nontrivial reorganization of terms, one obtains an inho mog eneous equation for J m (q 2 ), and 



a homogeneous one for m 2 (q 2 ), of the general form 29|, |30] 

J m {q 2 )[l + G{q 2 )} = 1 + tf JA^K^q,™ 2 ,^), 

/AdV. 
^-p/C 2 (M,m 2 ,A m ), (2.8) 



where d = 4 — e is the space-time dimension and p the 't Hooft mass. Note that the 
corresponding kernels are such that, as q — > 0, JCi^ik, q, m 2 , A m ) ^ 0; in fact, in Euclidean 
space, the solution of these equations furnishes J m (q 2 ) > and m 2 (q 2 ) > 0, as expected on 
physical grounds. 

The appearance of the factor [1 + G(q 2 )} on the lhs of Eq. (12. 8p stems from the fact 
that the corresponding equations are not derived from the standard SDE for the gluon 
propagator, but rather from its PT-BFM version; the advantages of this particular approach 
have been explained in detail in a series of articles [2ll-l23| . In the PT-BFM formalism the 
natural separation of the gluonic field into a "quantum" (Q) and a "background" (B) gives 
rise to an extended set of Feynman rules, and leads to an increase in the type of possible 



Green's functions that one may consider. In the case of the gluonic two-point function, in 
addition to the conventional QQ gluon propagator, A, two additional quantities appear: 
the QB propagator, A, mixing one quantum gluon with one background gluon, and the BB 
propagator, A, with two background gluon legs. It turns out that these three propagators 



are related by the all-order identities 24 



A(^) = [1 + G(q 2 )}A(q 2 ) = [1 + G(q 2 )] 2 A(q 2 ), (2-9) 

usually referred to as Background-Quantum identities (BQIs). 

The function G(q 2 ) is defined as the g^ v component of the special two-point function 

KM) = -i9 2 ^C A J ^A;(k)D(q-k)H ua (-q,q-k,k) 

= 9^G(1 2 ) + qj ^L(q% (2-10) 

q A 



where D ab (q 2 ) = 5 ab D(q 2 ) is the ghost propagator, and H ua is the gluon-ghost kernel [371.138]. 
Pertubatively, at one-loop, we have (Landau gauge, MOM scheme) 

l + G(g 2 ) = l + -^Clnf4V (2-11) 

Just as the usual quantum gluon propagator, the BB propagator A is also infrared finite, 
and must be parametrized in complete analogy with Eq. (12. 6p . namely 

A m 1 (g 2 ) = g 2 J m (g 2 )-m 2 (g 2 ), (2.12) 

and an exactly analogous formula holds for A (not used here). Then, one can establish that 
the BQIs hold individually for the kinetic and mass terms, i.e., 

J m (q 2 ) = [l + G(q 2 )] 2 J m (q 2 ), (2-13) 

m 2 (q 2 ) = [1 + G(g 2 )] 2 m 2 (g 2 ) (2.14) 

Use of Eq. (12.131) . together with Eq. ( 12. lip and Eq. ( 12. 5ft . furnishes the one-loop perturbative 
expression for J{q 2 ), namely 

J(g 2 ) = l + 6/ln^, (2.15) 

where b = 11CU/487T 2 is the first coefficient of the QCD /3-function. Evidently, as is well- 



known 



20| . the propagator A(q 2 ) absorbs all the RG logarithms, exactly as happens in QED 



with the photon self-energy 



III. EFFECTIVE CHARGE AND RG-INVARIANT GLUON MASS 

Let us next recall that, due to the Abelian Wis satisfied by the PT-BFM Green's func- 
tions, the renormalization constants of the gauge-coupling and of A(g 2 ), defined as 

A(gV) = ^V)Ao(<7 2 ), (3.1) 

where the "0" subscript indicates bare quantities, satisfy the QED-like relation 

Z g = Z A 1/2 . (3.2) 

As a result, the product 

d (q 2 ) ee glMq 2 ) = g 2 A(q 2 ) = d(q 2 ), (3.3) 

forms a RG invariant (/i-independent) quantity. As has been explained in the recent litera- 
ture 28|], d (q 2 ) may be cast in the form 



til" 

q 2 + m 2 (q 2 ) 



did') = „,,Z,n2V ( 3 - 4 ) 



with 

t(q 2 ) = 9 2 J m 1 (q 2 )i (3-5) 

and 

m 2 {q 2 ) = fh 2 {q 2 )J m \q 2 ). (3.6) 

Note that the two quantities defined above are individually RG invariant . The usual effective 
charge, a(q 2 ), is obtained from Eq. (13.5J) simply as a(q 2 ) = ~g 2 (q 2 )/ATt. At one-loop, 

_ g 2 1 

9 ( Q ^ = 1 + bg 2 In (q 2 /fi 2 ) = 6 In (q 2 /A 2 ) ' ^ 

where A denotes an RG invariant mass scale of a few hundred MeV. 

It is interesting to observe the analogy between the RG invariant mass defined in Eq. (J3.6P 
and the corresponding constituent quark mass, familiar from a plethora of studies on chiral 
symmetry breaking. Specifically, the quark propagator is usually cast in the form 

S-\p) = A{p 2 )i,-B{p 2 )l 

= A(p 2 )ty-M(p 2 n, (3.8) 



where I is the identity matrix, and the term A~ 1 (p 2 ) is often referred to in the literature 
as the "quark wave function". Note that, the quark acquires a dynamical mass (signaling 
the breaking of chiral symmetry) provided that B(p 2 ) is different from zero, and that the 
constituent quark mass, Ai(p 2 ) = B(p 2 )/A(p 2 ), is RG invariant . Then, it is natural to 
propose an analogy between the gluon and quark propagators; clearly, J m (q 2 ) corresponds 
to A(p 2 ), while m 2 plays exactly the role of B(p 2 ). Then, the division by J m (q 2 ) and A(p 2 ) 
gives rise, in both cases, to RG invariant masses, suggesting a close correspondence between 
m 2 (q 2 ) and M(p 2 ). 

Next, using the BQIs ( 12. 13ft and ( 12. 14ft to relate the components of A m (q 2 ) to the corre- 
sponding ones of A m (g 2 ), we get 

mV)£V) = ™V)-CV), (3-9) 

which finally furnishes a set of relations equivalent to (13. 5 p and (13. 6p , 

f(q 2 ) = g 2 [l + G(q 2 )}- 2 J m \q 2 ), (3.10) 

m 2 (q 2 ) = m 2 {q 2 )[l + G{q 2 )]- 2 J m l {q 2 ). (3.11) 

The basic Eqs. (13.51) and ( 13.61) . or alternatively Eqs. ( I3.10p and ( 13.111) . express the funda- 
mental quantities g 2 (q 2 ) and m 2 (q 2 ) in terms of J m (q 2 ), or J m (q 2 ) and G(q 2 ), respectively. It 
is therefore important to review the state-of-the-art in the determination of these quantities, 
both in the continuum as well as on the lattice. 

To that end, let us first focus on the quantities J m (q 2 ) and G(q 2 ), defined in the context 
of the conventional covariant gauges; in fact, we will specialize the discussion in the Landau 
gauge, where lattice simulations are usually performed. 

Let us say from the outset that lattice simulations of the gluon propagator A m (q 2 ) cannot 
furnish J m (q 2 ) without any additional input, for the simple reason that they provide the 
entire combination of J m (q 2 ) and m 2 (q 2 ), as appears in Eq. f |2.6|) . but not the individual 
components comprising it. The SDEs, on the other hand, when appropriately reorganized, 
give rise to two coupled integral equations, one for J m (q 2 ) and one for m 2 (q 2 ), as shown 
schematically in Eq. (I2.8p . In principle these two equations should furnish the exact behavior 
of J m {q 2 ) and m 2 (q 2 ); in practice, one is limited by the fact that the closed form of the kernels 
/Ci and JC 2 are only approximately known. In fact, /C 2 is better known than /C 1; mainly 
because the only unknown quantity that enters in /C 2 is the full three-gluon vertex T aiiu , 



whereas in /Q appears, in addition, the full four-gluon vertex. Even though neither of these 
two vertices is known, the STIs that the three-gluon vertex satisfies [see Eq. (14. 5)) ] furnish 
nontrivial information on the structure of its longitudinal components, thus providing certain 
field-theoretically motivated approximations for /C2; on the contrary, no such construction 
exists for the four-gluon vertex, which remains practically unexplored. 

The situation described above has motivated the advent of an approach which combines 
the information from both the lattice and the SDEs. In particular, instead of solving the 
system of Eq. (12. 8p . only the mass equation (involving the better known /C2) is solved, using 
as input in it the lattice data for A m (q 2 ). Then, once a solution for m 2 (q 2 ) has been obtained, 
one can use Eq. (12. 6 p to extract the approximate form of J m (q 2 ). 

Regarding the function G(q 2 ), notice that, quite remarkably, in the Landau gauge it 
coincides with the so-called Kugo-Ojima function [26|, |27J. In addition, there is a deep 
connection between the form factors G(q 2 ) and L(q 2 ) and the ghost dressing function 

F(q 2 )=q 2 D(q 2 ), (3.12) 

expressed through the identity [26|, [27J 

l + G{q 2 ) + L{q 2 ) = F-\q 2 ). (3.13) 



Since in d = 4 it is known that L(q 2 ) <C G(q 2 ) over the entire momentum range 27J, one 
finally arrives at the result 

l + GfcftwF-V), (3.14) 

which becomes an exact relation at q 2 = 0, since in this case L(0) = 0. 

From the above discussion it is clear that one cannot obtain direct information on J m (q 2 ) 
from the gluonic two point sector. Turning to the three-point functions, one might think 
that, given that J m (q 2 ) enters into the longitudinal c omp onents of the conventional three- 
gluon vertex, as first demonstrated by Ball and Chiu J35[, a judicious combination of them 
could project it out. However, the problem in this case is the "contamination" from the 
various form- factors comprising H^ u , which enter nontrivially in the various expressions. 

It turns out that the PT-BFM three-gluon vertex lends itself for this type of analysis. 
Given that it satisfies QED-like Wis instead of STIs, there is no ghost sector contributions: 
the longitudinal form factors can be expressed in terms of the J m only. To be sure, the 
transverse form factors enter in general, and are unknown. But, as we will see in Sec. |V] there 

10 



is a simple kinematic limit, studied usually in the lattice simulations of the (conventional) 
vertex, which eliminates completely the transverse components, and projects out solely the 
first derivative of q 2 J m (q 2 ). It is therefore clear that a possible lattice extraction of J m (g 2 ) 
would provide an immediate determination of ~g 2 {q 2 ) from Eq. ( 13. 5p . 

On the other hand, regarding the prospects of extracting the RG invariant gluon mass 
from Eq. (13. 6p or ( 13. lip , one may envisage two basic scenarios [again assuming independent 
knowledge of J m (<l 2 )]'- 

(i) A m (q 2 ) has been simulated on the lattice. Then, it is direct to extract m 2 (q 2 ) from 
Eq. (I2.12p . since (Euclidean space) 

m 2 (q 2 ) = A m \q 2 )-q 2 J m (q 2 ), (3.15) 

or, using Eq. ( I3.6p . 

m 2 (q 2 ) = A- m \q 2 )J m \q 2 )-q 2 . (3.16) 

(ii) A m (q 2 ) has not been simulated on the lattice. Then, to proceed further, one must 
necessarily employ the BQIs; this is rather feasible, because, as mentioned earlier, 
the function G(q 2 ) has been simulated on the lattice. Specifically, one may use the 
conventional A m (q 2 ) obtained from the lattice together with G(q 2 ) to build A m (q 2 ) 
by means of the BQI. Then one may return to the previous case (i), and substitute it 
into Eqs. (13151) and <KW\ . 

This general discussion motivates a systematic study of the PT-BFM three-gluon vertex. 
To that end, in the next section we will present a brief reminder on the structure and 
general properties of the B 3 vertex, while in Sec. |V] after summarizing the current prospects 
of simulating nonperturbative PT-BFM Green's functions on the lattice, we will consider 
how J(q 2 ) maybe extracted from a possible lattice simulation of this vertex. 

IV. THE PT-BFM THREE-GLUON VERTEX 



The gauge- invariant three-gluon vertex T afJ-l/ has been first considered in 15| . where its 
one-loop construction was carried out by means of the PT, and its basic WI was derived 

n n 

(see also [19(). It was further studied in [3Jj, with particular emphasis on the special rela- 



tions between gluonic, fermionic, and scalar loop contributions, and has been revisited very 
recently in [39(, using string-inspired techniques. 

11 





FIG. 1: The one-particle irreducible B 3 vertex and the SDE it satisfies. Background legs are 
indicated by the small circle at their end; black circles correspond to 1PI Green's functions, white 
circles represent connected functions, while gray circles indicate SDE kernels. The corresponding 



symmetry factors may be found in 



M- 



A. General properties 



The equivalence between PT (or the "generalized PT") and BFM allows for a concise 
field-theoretic definition of this vertex, as vacuum expectation value of the time-ordered 
product of three background gluons. In particular, denoting the full connected three-point 
function by (momentum space) 



nabc 



(q,r,p) 







T[A*MA b Jr)AZ(p)] 







(0|0) 



one defines 



Q%,{<l,r,p) = A^(g)A»,(r)A£(p)r^(?,r,P) 



(4.1) 



(4.2) 



where T alJiV (q,r,p) is the amputated three-point function. 

Exactly as happens with the conventional three-gluon vertex T 



aiil/ , the T aill/ is fully Bose 



symmetric. This is to be contrasted with the BQ 2 vertex, usually denoted by T aiiu , which 
is Bose symmetric only under the interchange of its two quantum legs. The SDE satisfied 
by T aiiu is shown in Fig. CD Note that the tree-level expressions for T afJLV and T afll , coincide: 



r 2L(?> r >P) = ri°L(g, r,p) = (q- r) v g ail + (r - p) a g lll/ + (p - q)^g av . 



(4.3) 
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For the purposes of the present work, the most important property of T afll/ is the fact 
that it satisfies the completely Bose-symmetric set of Abelian Wis 

q a T afl „(q,r,p) = p 2 J(p 2 )P^(p) - r 2 J(r 2 )P„„(r), 



rT 



aiiv\Hi 



r,p) = q 2 J(q 2 )P a v(q)-p 2 J(p 2 )P a „(p), 



p»r a , u (q,r,p) = r 2 J{r 2 )P a ,{r) - g 2 J(g 2 )P Q/i (g). 



These simple Wis are to be contrasted with the STIs satisfied by T atiu , namely 



(4.4) 



q a T a ^(q,r,p)=F(q 2 )[p 2 J(p 2 )P«(p)H a ,(p,q,r)-r 2 J(r 2 )P«(r)H ai/ (r,q,p)] (4.5) 



and cyclic permutations. The tensorial decomposition of H Vfl is given by [35 ] 



■Lli/jj,\p, Tf q) QiivQ'qrp i/j,Qu"qrp T qfiPuCqrp "T qvPfi.O'qrp "T Pp!Pv& 

where a qrp is short-hand notation for a(q,r,p), etc. 

An immediate consequence of Eq. (14. 4 p is the QED-like relation 

Z\ = Za- 



qrp ■ 



(4.6) 



(4.7) 



between the wave-function renormalization for A, introduced in Eq. f 13 .11) . and the vertex 
renormalization defined as 

Z 1 f^(q,r,p) = f a R ^(q,r,p). (4.8) 

In addition, one may extract from the set of Wis given in Eq. ( 14. 4 p the expression of 
the vertex T ativ in the kinematical limit r — > 0, which will be the relevant momentum 



configuration in the subsequent analysis (for a related analysis, see also 40]). To show that, 
consider the Taylor expansion of a function f(q, r,p) around r = (and p = —q). In general 
we have 



f(q,r,p) = f(q,0,-q)+r^ 



_d_ 



f(q,r,p)} +0(r 2 



(4.9) 



r=0 



where the Lorentz structure of the function / has been suppressed. Specializing this result 
to the second WI in Eq. (14.4[) . one finds 



r»r a Uq,r,p) = rT«™(g,0, -q) + 0{r 



d 



(q + r) 2 J(q + r)P au (q + r) 



+ 0{r 2 



(4.10) 



r=0 
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where we have use the fact that the zero-th order term of Eq. (14. 9 j) vanishes in this case. 
Thus, equating the coefficients of the terms linear in r M , one obtains the relation 

f d r 
IV(g, 0, -g) = -l-^[(q + r) 2 J(q + r)P a ,(q + r) 

d 



r=0 
2 ?/„ 2\ 



q 2 J(q')P av (q) • (4.11) 



<9g^ 

Interestingly enough, the above relation is formally equivalent to the well-known QED text- 
book result 



d_ 

dp 
obtained from the WI 



T,(0,-p,p) = -—S- 1 (p), (4.12) 



<7T M (<7,-g-p,p) = S-\p) - S-\q + p), (4.13) 

relating the photon-electron vertex with the electron propagator. 

Next, the derivative of Eq. ( 14. lip can be easily evaluated by applying the formula 

O Q Q \ 

;Pau{q) = 2q»^f ~ —{g^aqv + g^,q a ) 



= --\q a P^{q) + quP^a(q)}, (4.14) 

yielding the final result 

r aM v(g, 0, -q) = J{q 2 ){g aii qy + gu^q a - 2g au q^) - 2q 2 J\q 2 )q fl P au (q) 

= J(q 2 f { *lM^,-q)-2q 2 J\q 2 )q,P av (q), (4.15) 

where the prime indicates the derivative with respect to q 2 . Note that, due to Bose symmetry, 
the procedure described above can be applied exactly in the same way for the remaining 
Wis, in order to obtain the kinematical limits q,p — > of the vertex T afJiU . Also, we observe 
that the tensorial structure of the full vertex T atiu in this particular kinematical limit is not 
exhausted by the term proportional to the tree level expression. 

B. The pole part of the three-gluon vertex 

As has been explained in a series of recent works, a crucial condition for obtaining an 
infrared finite gluon propagator, without interfering with the gauge invariance (or the BRST 
symmetry) of the theory, is the existence of a set of special vertices, to be generically denoted 
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by V, that are completely longitudinal and contain massless poles. The dynamical origin 
of the aforementioned poles is purely nonperturbative: for sufficiently strong binding, the 
mass of certain (colored) bound states may be reduced to zero [4lN45] . 

The role of these vertices is two-fold. On the one hand, thanks to the massless poles 
they contain, they make possible the emergence of an infrared finite solution out of the SDE 
governing the gluon propagator; thus, one invokes essentially a non-Abelian realization of 
the well-known Schwinger mechanism 46|,|47j]. On the other hand, these same poles act like 
composite Nambu-Goldstone excitations, preserving the form of the STIs of the theory in 
the presence of a gluon mass. 

Specifically, in order for the Wis to maintain the same form before and after mass gen- 
eration, the effective substitution 



a- V) = q 2 J(q 2 ) — ► a,; 



-1( J2\ 



q 2 J m (q 2 )-fh 2 (q 2 ), 



(4.16) 



implemented by the mass generation at the level the gluon propagator, must be accompanied 
by the simultaneous replacement of the vertex 29] 



r = r m + v. 



(4.17) 



Then, since 

q a F ma nv(q,r,p) = p 2 Jm{p 2 )P^u{p)-r 2 J m {r 2 )P^u{r)., 
q a V ai , u (q,r,p) = m 2 (r 2 )P^(r) - m 2 (p 2 )P^ v (p), 

one finds that the corresponding WI satisfied by r" would read 



9°T' (q,r,p) 



Fm{q,r,p)+V{q,r,p) 



afiv 



[p 2 J m (p) ~ m 2 {p 2 )]P» v {p) ~ [r 2 J m {r 2 ) - m 2 {r 2 )]P^{r) 



A™V)^(p)-A m V)iV(r), 



(4.18) 
(4.19) 



(4.20) 



which is indeed the first of the identities in Eq. (14 .4|) . with the aforementioned replacement 
A -1 — > A^ 1 enforced. 

The closed expression for V afliy may be reconstructed from the Wis it satisfies (namely 
Eq. (14.191) and its cyclic permutations), together with the condition of complete longitudi- 
nality, i.e., 

P a ' a (q)P»'»(r)P»'»(p)V a ,,, u ,(q,r,p) = 0; (4.21) 
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it reads, 

V a , u (q,r,p) = ^[m 2 (r 2 ) - m 2 (p 2 )]P^r)P pu (p) + c.p., (4.22) 

where c.p. denotes "cyclic permutations". As we will see in the next section, due to its 
completely longidudinal nature [viz. Eq. (J4.21J) ]. the contribution of this pole vertex to the 
typical lattice ratio, given in Eq. f 15 . 3 1) . vanishes in the Landau gauge. 



C. The longitudinal form factors 

The complete closed form of T m is not known; its longitudinal part, however, may be 



reconstructed from the Wis that T m satisfies, following rather standard procedures 35 ] . 
Specifically, one begins by separating the vertex into the "longitudinal" and the (totally) 
"transverse" parts, 

fTfer.p) = f^wfor.p) +C5(?,r,p), (4.23) 

where the component T m ^ satisfies the WI of Eq. (14.181) (and its permutations), whereas 
^C^(g,r,p) = r^ t) {q,r,p)= Pv Y a ^ t) {q,r,p) = 0. 

The longitudinal part is then decomposed into 10 form factors X^, according to 

10 

f a m ^(q,r,p) = J2MQ,r,p)£r, (4-24) 



with the explicit form of the tensors t given by [31| 

t^ v = (q - r) u g a ^ t^ = -p v 'g a » £^ v = (q - r) u [q^r a - (q ■ r)g a ^] 

t/ v = (r - p) a g^ Q" v = -q a g^ 1%^ = (r - p) a [r v p» - (r • p)g^} 

t^ = {p- qYg av t/ v = -r»g av t^ v = (p - qY[p a q u - (p • q)g au ] 

^iq U = q u r a p fl + qf^^p . 

Then, the WI of Eq. (14. 18ft and its permutations give rise to an algebraic system for the Xj, 
whose solution reads, 

X 1 = ^[J m (? 2 ) + ^(r 2 )], X 2 = i[J m (g 2 )-J m (r 2 )], x 3 - J ™^ ~ J ^ 

X, = ^[J m (r 2 ) + J m (P% X b = l -[J m {r 2 )-J m {p 2 )l X 6 

X 7 = \[J m {p 2 ) + J m {q% X 8 = ^[J m (p 2 )-Jm(q% X 9 

X 10 = 0. (4.26) 
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q 2 - 

•J my ) 


~ Jm{p 2 ) 


Jm(p 2 ) 


- p 2 

- Jm(q 2 ) 


P 2 


-q 2 



Thus, the longitudinal form factors of r^Jg,r, p) involve only the quantity J m , since it is 
the only quantity that enters on the rhs of Eq. ( I4.18p . Instead, the correspon ding expressions 
for the form factors of the conventional (Q 3 ) vertex, Y afiu , first derived in [35j . contain, in 
addition, the ghost dressing function F and the various form-factors comprising the gluon- 
ghost kernel H^ u . For example, the closed form of X 7 (to be employed in the next section) 
is given by 

X 7 (q, r,p) = -{2[F(q)J m (p)a rqp + F(p)J m (q)a rpq ] + r 2 [F(p)J rn (r)b qpr + F (q) J m (r)b pgr ] 
+ (<? 2 - P 2 )[F{r)J m {q)b pr q + F(q)J m (p)b rqp - F(r)J m (p)b qrp - F(p)J m (q)b rpq ] 
+ 2(qr)F(p)J m (q)d rpq + 2(rp)F(q)J(p)d rqp }. (4.27) 

Note in addition, that, unlike X 10 , the corresponding X w does not vanish. 

Finally, the (undetermined) transverse part of the vertex is described by the remaining 
4 form factors K, 



■ i-i 



f a m % ) (q,r,p) = J2U^r,p)tr, (4.28) 



with the completely transverse tensors t l given by 

t^ v = [(q ■ r)g a » - q»r a ] [(r • p)q u - (q ■ p)r v ] 
t%" = [(r ■ p)g^ - r V] [(p ■ q)r a - (r • q)p a ] 

tT = i(p-q)g ua -p a q u ]{(q-r)p»-(r-p) q »] 

tf v = g^[(p ■ q)r a - (r • q)p a ] + g afl [(r ■ p)q v - (q ■ p)r v ] + g au [{r ■ q)p fl - (r • p)tf] 

+ p a q^r u - r a p^q u . (4.29) 

It turns out that in the limit r — > all the transverse tensors in Eq. (I4.29P are zero, so 
the transverse part of the vertex vanishes in this limit. On the other hand, for the same 
limit, only the following longitudinal tensors given in Eq. ( I4.25P survive 



(4.30) 



t^ v = q"g a » t^ v = q u g Q » tf v = q a g^ 

t K»v = _ q(Xg „ v t ^v = _ 2q n g av %*» = -2q 2 q^P a %q), 

with the associated form factors 

Xi = \[J{q 2 ) + J(0)], X 2 = i[J(g 2 )-J(0)], X 4 = i[J(0) + J(g 2 )], 

X, = \[J{0) - J(q% X 7 = J(q 2 ), X 9 = J\q 2 ), (4.31) 
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corresponding to the limit r — > of those appearing in Eq. (I4.26p . Thus, using these results, 
it is elementary to show that one is able to reproduce the expression given in Eq. (I4.15p . 

V. LATTICE PROSPECTS 

In this section we first discuss the theoretical possibilities of simulating PT-BFM Green's 
functions on the lattice. Then, we consider the relevant lattice quantity, and derive its 
general expression in terms of the various (longitudinal and transverse) form factors. Next, 
we show how the effective charge and gluon mass may be reconstructed from the data 
obtained in a standard kinematical limit. In addition, we present a numerical study on the 
propagation of the (modeled) data errors into the effective charge and gluon mass. Finally, 
exploiting some basic field-theoretic properties, we relate the conventional and BFM lattice 
quantity at the origin. 

A. BFM on the lattice 



Within perturbation theory, the BFM was formulated to all orders on the lattice in 48] 



see also the early work of Gross and Dashen 49J); however, its nonperturbative implementa 



tion has been pending for quite some time. A possible nonperturbative formulation, which 
evades the well-known "Neuberger 0/0 problem" [50J, has been only recently introduced 



in 



321 ]. through a reformulation of the BFM method in terms of canonical transformations 



so that dynamical ghosts are not needed). This led to the proposal of the gauge fixing 



functional 
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33] 



F[g} = -Jd i xTr(A^-A fl ) 2 , (5.1) 

which upon minimization on the group elements g provides the background Landau gauge 
condition V^(A 9 — A^) = 0. Then, on the minimum of this functional, the mapping given 

by 

A -> A 9{A ' A) -A = Q^ (5.2) 

defines the action of a canonical transformation on the gauge fields, constituting a non- 
perturbative generalization of the familiar splitting into a background and a quantum field 
5JJ. The inverse of this mapping amounts to a gauge transformation, which can be used 



to determine the quantum field Q corresponding to the gauge configuration minimizing the 
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gauge-fixing functional ( 15. ip . The relevant correlators (e.g., in the two- and three-point 
functions) may then be constructed in terms of this particular field, and lattice simulations 
of the corresponding discretized version can be then carried out. 

Preliminary numerical simulations for a variety of explicit backgrounds have been per- 
formed in 33(, indicating that the functional (15. ip furnishes a convergence rate comparable 
to that of the zero background case. Note, however, that in the PT-BFM formulation one 
considers the background field as an external yet unspecified source, to be set to zero after 
taking the appropriate derivatives of the vertex functional. Therefore, in order to properly 
simulate this procedure on the lattice, one ought to introduce a suitable dependence of the 
background field A on a parameter, which, upon variation, would smoothly turn it off. Work 
in this direction is already in progress; at the moment, we are not aware of any theoretical 
obstruction that would prevent the computation of PT-BFM Green's functions on the lattice 



by techniques similar to those discussed in [52]. 



B. The basic lattice quantity 

Let us therefore assume that the three-point function defined in Eq. (14. ip may be in- 
deed simulated on the lattice, following the procedure briefly outlined above. Then, as is 



customary, in the Landau gauge, one considers the following ratio 36] 



R(q,r,p) = — -, (5.3) 

V(q,r,p) 

with the numerator and denominator given by 

N(q,r,p) = T^l v (q,r,p)g puT (q,r,p), 

V(q,r,p) = r^(g,r,p)P^(g)P^(r)P-(p)r(°J r (g,r,p)A m (g)A m (r)A m (p), (5.4) 

where a common color factor cancels out in the ratio. Note the index "m" in the full 
propagators, indicating the nonperturbative generation of a gluon mass. In addition, and 
according to the discussion in the previous section, the substitution given in Eq. ( I4.17P 
must also be implemented. Then, inserting Eq. ( 14.21) into Af(q,r,p), together with the 
aforementioned substitution, we see that (i) the product A m (g)A m (r)A m (p) cancels out 
when forming the ratio of Eq. (15.31) . and (ii) any reference to the pole vertex V paT disappears, 
due to the longitudinality condition Eq. (I4.2ip that it satisfies. 
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Thus, the numerator and denominator become 

Af(q,r,p) = r^( g ,r,p)P;(g)P^(r)P;(p)r^(g,r,p), 

V{q,r,p) = T^ u (q,r,p)P a ''(q)P^(r)P^(p)T^ T (q,r,p). (5.5) 

When one decomposes the full three-gluon vertex into a longitudinal and a transverse 
part, as in Eq. (I4.23p . the numerator in Eq. (15.51) becomes 

Af(q,r,p) =N(e)(q,r,p)+J\f( t) (q,r,p), (5.6) 

with 

Af m (q,r,p) = T^(q,r,p)P^q)P^r)P:(p)f^ t) (q,r,p). (5.7) 

Then, the denominator is given by 



V(q, r,p) = 4 r P _ ( r W [ 3 ( g 2 r 2 + g2p2 + ^ + ^2 _ ^2^ ^ 



the longitudinal part of the numerator by 

N(e)(q,r,p) = 4 F 2 ^ * } |[3g 2 r 2 - {q-p){p-r)]A 1 + [3r 2 p 2 - {p-q){q-r)]A 2 

+ [3<? 2 p 2 - (<?-^)(^-p)]^3 + l(q- r )( r -p)(p-q) - q 2 r 2 p 2 ]A 4 \ , (5.9) 

and its transverse part by 

Af {t) (q,r,p) = 2[rV - [rpf] {[3(g-r) -p 2 ]^ + [3(r-p) - g 2 ]F 2 

+ [3(g-p) - r 2 ]? 3 + 6F 4 } • (5.10) 

Note also that the identity 

(g-r)(r-p) + (r-p)(p-q) + (p-q)(q-r) = q 2 r 2 — (g-r) 2 = g 2 p 2 — (q-p) 2 = r 2 p 2 — (r-p) 2 , (5.11) 

which can be easily proved using momentum conservation, has been employed in deriving 
the above expressions, and that we have introduced the notation 

A 1 = X 1 -(q- r)X 3 ; A 2 = X 4 - (r-p)X 6 ; 

A 3 = X 7 -(p-q)X 9 ; A 4 = -X 3 - X 6 - X 9 . (5.12) 

It is important to emphasize that the expressions given in Eqs. (15. 9p and (J5.10P carry 
over directly to the case of the conventional three-gluon vertex T^ , simply by converting 
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the hatted quantities to normal ones. Of course, as already mentioned, the functional form 
of the corresponding Xi comprising the Ai is significantly more complicated. Note also 
that even though X 10 does not vanish, it still does not contribute to Eq. (15. 9p . because 
P«(q)PZ(r)P»(p)l^ = 0. 

It turns out that the physical quantity of interest may be extracted from R by employing 
a very common kinematic choice. Specifically, to begin with, as is customary in lattice 
studies, we will describe the ratio R in terms of the modulo of two independent momenta 
(say, q 2 and r 2 ) and the angle formed between them; thus R = R(q 2 , r 2 , 0) . Then, the 
quantity of interest corresponds to the case R(q 2 , 0, vr/2), which is a special case of the 
so-called "orthogonal configuration", namely R(q 2 ,r 2 ,7r/2). 

In this latter configuration we have, 

p 2 = q 2 + r 2 ; q-r = 0; q-p=—q 2 ; r-p=—r 2 , (5.13) 

and, therefore, the relevant quantities reduce to 

V(q, r, vr/2) = ^^[3(g 4 + r 4 ) + 8gV], (5.14) 

and 

AfW(q, r, vr/2) = — ^-\2q 2 r 2 A x + 3r 2 (g 2 + r 2 )A 2 + 3q 2 (q 2 + r 2 )A 3 - q 2 r 2 (q 2 + r 2 )A 4 ], 
q z + r z 

A/" (t) (g,r,vr/2) = -2g 2 r 2 [(g 2 + r 2 )F 1 + (3r 2 + g 2 )F 2 + (3g 2 + r 2 )l > 3 + 6i > 4], (5.15) 

with the (suppressed) arguments of the form-factors A4 and Y { correspondingly adapted to 
the particular kinematic configuration chosen. 

At this point, if, in addition, we set r 2 = 0, then the transverse term vanishes, 
A/" (t) (g, 0, 7r/2) = 0, and we obtain 

R(q 2 , 0, vr/2) = A 3 (q, 0, tt/2) = X 7 + q 2 X 9 , (5.16) 

so that (we only indicate the q 2 in the argument of R) 

R(q 2 ) = [q 2 J m (q 2 )]', (5.17) 

where, as before, the prime indicates derivatives with respect to q 2 . Let us point out that this 
particular result may be derived directly from Eq. ( 15.3[) . by substituting in it the expression 
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for T afll/ given in Eq. (j4.15p . Specifically, 

R (q ,r, p) = J {f) - VTM { Ei^^fgmSW J . (5 , 8) 

In the orthogonal configuration (g-r = 0), we have that q a P£(r) = q**; then, using Eq. (I5.14p 
and setting r 2 = 0, it is easy to show that the quantity in curly brackets is equal to {— |}. 
Then, simple integration of Eq. (I5.17P yields 

q 2 Jm(Q 2 )= / dp 2 R(p 2 ) + C, (5.19) 

Jo 

and, assuming that both J m {q 2 ) and R{p 2 ) are finite for all values of the momentum, we see 
that the integration constant must vanish, C — 0. Thus, finally, one obtains the relation 

Jm{q 2 ) = ~ 2 [ q dp 2 R(p 2 ). (5.20) 

Let us next renormalize this result within the MOM scheme, denoting the final answer 
by Jm (q 2 ). If we impose the standard MOM condition Jm (/x 2 ) = 1, at some arbitrary 
momentum scale /x, then, we have that 

%V) = i^- (5-21) 

7 (a 2 ) 

However, let us point out that, strictly speaking, due to the presence of the gluon mass, 
this last normalization condition imposed on J m cannot be enforced simultaneously with 
the corresponding MOM condition for A~ x (g 2 ), namely A~ 1 (/i 2 ) = 1. Indeed, since, for any 
arbitrary /i, A~ 1 (/i 2 ) = fi 2 J m (fi 2 ) +m 2 (/x 2 ), if at a given \i we impose that J m (fi 2 ) = 1, then, 
automatically, at the same /i, A~ 1 (/i 2 ) = /i 2 [l + m 2 (/i 2 )//i 2 ]. Note however, that unless one 
chooses to push the value of /x very deep in the infrared, this discrepancy is numerically 
immaterial; for example, when one renormalizes A m at /x = 4.3 GeV, the value of the gluon 
mass at the origin is m(0) = 1 GeV, making the ratio m 2 (0)//x 2 of the order of 5%. Of course, 
this estimate is just an upper bound for the relevant ratio m 2 (/x 2 )//x 2 , which in reality is 
significantly smaller, since the function fh 2 (q 2 ) is decreasing rather rapidly (see, e.g., the 
inset of the left panel of Fig. [2]); in fact, at /x = 4.3 GeV, the gluon mass is practically 
negligible. 
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FIG. 2: The expected shape of the ratio R(q 2 ) (right panel), obtained starting from the quenched 
lattice data for the SU(3) gluon propagator A m (q 2 ) (left panel), and the corresponding solution 
of the mass equation m 2 (q 2 ) (left panel, inset). 

C. Modelling the error propagation 

In order to understand how a possible lattice signal for the quantity R(q 2 ) of Eq. (I5.17P 
may provide direct information on the effective charge a(q 2 ) (and the possible caveats asso- 
ciated with such a determination), we perform in what follows a detailed numerical study. 
Specifically, starting from the knowledge of the conventional quenched lattice propagator 
A m (q 2 ) and ghost dres sing function F(q 2 ) [3J, together with the associated solution of the 
mass equation m 2 (q 2 ) [30], one can reconstruct first J(q 2 ), and next obtain the expected 
shape of R(q 2 ) through 



R{q 2 



A m V)-^V) 
F 2 (q 2 ) 



(5.22) 



where we have used Eq. f)2.6p in Euclidean space. 

The resulting curve (shown in Fig. |2]), which will be referred to as the 'expected' result, 
can be parametrized to a high precision by the function 



R{q 2 ) = A 2 + 



A t -A 2 



with best fit parameters corresponding to the values 



A± = 0.083; A 2 



5.150; 



% 



7.156 GeV; 



x 



0.836. 



(5.23) 



(5.24) 



The above functional form of the expected behavior of R{q 2 ) is rather useful, because 
it allows for a systematic analysis of how uncertainties, simulated through deviations of 
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FIG. 3: (color online) Numerical study of how uncertainties in the determination of R(q 2 ), 
parametrized according to Eq. H5.23H . reflect into the determination of J, J as well as the effec- 
tive charge a(q 2 ). In particular, we compare the 'expected' results with the results reconstructed 
from the parametrization (|5.23p for the following parameter values: (first row) the best fit param- 
eters (|5.24p ; (second row) A\ fixed at the expected value R(0) = 0.0562, and, for the remaining 
coefficients the refitted values A2 = 5.246, qo = 7.279 and x = 0.8167; (third row) Ai and x fixed 
at the values R(0) and 0.92 respectively, and A2 = 4.632, qo = 6.234 GeV ; (fourth row) A\ and x 
fixed at the values R(0) and x = 0.72 respectively, and, finally, A2 = 6.251 qo = 9.21 GeV. 

the fitting parameters from their "optimal" values, can influence the reconstruction of the 



24 





0.14 




- 




"""""j^ m 2 reconstructed 




0.12 




- 


CM 

> 

CD 


0.1 
0.08 


'•\ 


- 


,. — * 








o> 


0.06 




- 


£ 


0.04 




- 




0.02 







- 




\ 


/VX/^/VV, 








*/-/v!ffi, 











0.001 



0.01 



0.1 

q 1 [GeV 2 ] 



10 



FIG. 4: The reconstructed running mass starting from R(q 2 ) as given in Eq. (I5.23h when the param- 
eters chosen are the best fit ones of Eq. (|5.24p . Notice the problems that afflict the determination 
of the UV mass tail (shaded region). 

effective charge (and gluon mass, see below). 

The results of this study are shown in Fig. [3j As can be seen, the most sizable devia- 
tion between the 'expected' and the reconstructed results occurs when uncertainties in the 
determination of -R(O) are sizable (Fig. [31 first row). This is mainly due to the fact that 
the 'expected' value for R(0) turns out to be of C(10~ 2 ), and moderate deviations imply a 
considerable effect in the determination of J m through Eq. (J5.20P ; this, in turn, translates 
into a large variation of the effective charge, given that a ~ J" 1 , see Eq. (13. 5p . 

Specifically, in the first row of Fig. |3] we show the results for R, J, J and a obtained 
starting from Eq. (I5.23P with the best fit parameters ( I5.24P ; since in the parametriza- 
tion (I5.23p . R(0) = A 1 , such fit overestimates the value at the origin; this error propagates 
in the determination of a reconstructed effective charge, which comes out suppressed in the 
IR with respect to the 'expected' value. 

The following three rows in Fig. [3] show the effect of uncertainties in the determination 
of the exponent x (which controls the overall shape of the R curve), once the behavior at 
the origin has been fixed to its 'expected' value, that is A\ is fixed to the value R(0). As 
can be seen, the effect is significantly milder, and one can always reconstruct the effective 
charge to a high degree of accuracy. 

Finally, let us focus our attention to the mass m 2 (g 2 ); as we will see, its reconstruction 
from this particular type of lattice measurements is especially subtle. The main difficulty 
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FIG. 5: The ghost-gluon vertex and the Taylor kinematics. 

stems from the fact that the extraction of m 2 (q 2 ) proceeds by means of the relation 

mV) = A m V)-g 2 J m (g 2 ). (5.25) 

Now, whereas both terms on the rhs increase in the UV, on theoretical grounds we know that 
m 2 (q 2 ) must decrease rather rapidly; this, in turn, implies that a rather delicate cancellation 
between these aforementioned two terms must take place. This cancellation, however, is very 
likely to be distorted by the reconstruction procedure, especially for high values of q 2 . This 
particular problem is shown in Fig. HJ where we plot the reconstructed mass for R(q 2 ) 
provided by Eqs. f)5.23|) and (15.24)) . As can be seen there, while the determination in the 
IR is reasonably accurate, the tail is seriously distorted, displaying even negative regions. 
This characteristic pathology persists (with various degrees of intensity) in all cases analyzed 
in Fig. [3J 

D. Relating R(0) with R(0) 

Let us now consider the conventional R{q 2 ), obtained from Eq. ( 15.16p through the direct 
substitution X 7j g — > X 7) g. Obviously, the complicated structure of X 7 and Xg [see, e.g., 
Eq. (14.27)) ] . infested by the unknown form factors of the ghost-gluon kernel, makes their use 
for arbitrary q 2 impractical. However, when q 2 = the corresponding expressions simplify 
substantially, providing a fairly simple expression for R(0). 

Specifically, after implementing X79 — > A^g in Eq. (15.160 . let us set q 2 = 0, to obtain 

J R(0) = F(0)J m (0)a(0,0,0), (5.26) 

where Eq. ( 14.271) has been employed. 

At this point one may invoke Taylor's theorem 53( in order to determine, under mild 



assumptions, the value of a(0, 0, 0). To that end, consider the general Lorentz decomposition 
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for the ghost-gluon vertex, 

r M (r,g,p) = B 1 (r,q,p)q fl - B 2 (r,q,p)p fl . (5.27) 

Notice that, setting B{ = and B 2 = 1, one recovers the tree- level value of the ghost- 
gluon vertex T^ = —p p . In the Taylor kinematic configuration, see Fig. |5l Eq. ( I5.27P 
becomes 

r„(0, -p,p) = -[Bi(0, -p,p) + S 2 (0, -p,p)] Pll . (5.28) 

Taylor's theorem states that 

B 1 (Q,-p,p) + B 2 (Q,-p,p) = l, (5.29) 

to all-orders in perturbation theory. 

On the other hand, it is well-known that T /i (r, q,p) can be obtained from the contraction 

r M (r,g,p) = -p u H Pll (p,r,q) 

= -[a qrp + (q-p)b qrp + {q-p)d qrp + p 2 e qrp ]p fM - [{q-p)b qrp + p 2 c qrp ]q^, (5.30) 

which, in the Taylor kinematics reduces to 

F»(0,-p,p) = -{a(-p,0,p) -p 2 [c(-p,0,p) + d(-p,0,p) -e(-p,0,p)]}p M . (5.31) 

Thus, equating Eq. (I5.28J) with Eq. ( 15.311) and using Eq. (15.291) . we deduce the constraint 

a(-p,0,p) -p 2 [c(-p,0,p) +d(-p,0,p) -e(-p,0,p)] = 1. (5.32) 

Finally, taking the limit p — > in Eq. (I5.32p . and assuming that c,d, and e are regular 
functions in that limit one obtains that a(0, 0, 0) assumes its tree-level value, 

a(0,0,0) = l. (5.33) 

Consequently, Eq. f l5.26f) becomes 

R(0) = F(0) J m (0), (5.34) 

whereas, in the same limit, Eq. (15.1 7p gives directly 

R(0) = J m (0). (5.35) 
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It is now relatively straightforward to derive a simple relation relating -R(O) with R(0). 
Specifically, use of Eq. ( 12. 13[) and the relation (I3.13P at q 2 = 0, yields 

J m (0) = [1 + G(0)] 2 J m (0) = F- 2 (0) J m (0), (5.36) 

so that finally 

5^ = F 3 (0). (5.37) 

R(0) 

Since F(0) > 1, the determination of -R(O) from -R(O) might be error prone; in the case 

studied in the previous section, at /i = 4.3 GeV one has F(0) = 2.86, which furnishes 

-R(O) = 1.58 for the initial value -R(O) = 0.067. Therefore, given also the sensitivity of 

the effective charge to the value of i?(0), one should ideally proceed the other way round, 

limiting the possible values of -R(O) through an independent measurement of R(0). 

VI. CONCLUSIONS 

In this work we have explored the possibility of determining the complete momentum 
evolution of the QCD effective charge from a special kinematic limit of the three-gluon 
vertex corresponding to three background gluons. Given that within the BFM quantization 
scheme the (background) gauge invariance is preserved, the aforementioned vertex satisfies 
linear Wis, with no reference to the ghost sector; within the PT the same property emerges 
naturally, after the systematic rearrangement of an appropriate observable, following the 
standard pinching rules. Consequently, and in contradistinction to what happens in the 
case of the conventional three-gluon vertex, the longitudinal form factors of this PT-BFM 
vertex may be expressed exclusively in terms of the background gluon wave function J m . By 
virtue of the dynamically generated gluon mass, this latter quantity, as well as the physical 
effective charge defined from it, are infrared finite and free of any divergences related to the 
perturbative Landau pole. 

Particularly interesting in this context is the possibility of simulating the (Landau gauge) 
PT-BFM vertex on the lattice. To that end, after briefly reviewing the general theoretical 
feasibility of such a task, we have focused on the specifics of the vertex simulation, with 
special emphasis on the relevant kinematic limit that projects out the desired quantity 
J m (q 2 ). In addition, a preliminary numerical analysis suggests that the extraction of the 
effective charge is relatively insensitive to the numerical uncertainties that may infest an 
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actual lattice simulation. On the contrary, the reconstruction of the running gluon mass 
turned out to be considerably more subtle; specifically, instead of being positive-definite and 
monotonically decreasing, the corresponding curve displays unphysical fluctuations past 1 
GeV, due to the distortion of delicate numerical cancellations. 

The results presented in this article are expected to contribute to the collective effort 
dedicated to the deeper comprehension of the nature and properties of the QCD effective 
charge. In fact, the possibility of probing directly the value of u(q 2 ) at the origin, through 
the corresponding lattice extraction of J m (0), is intriguing, and may serve as testing ground 
for various alternative pictures. In particular, the value of a(0) is directly related to the 
notion of the "QCD conformal window", appearing in studies based on the AdS/QCD 



correspondence [13] . 



It is important to emphasize that the nature of the observable R, coupled with the 
fact that we work in the Landau gauge, results in the total annihilation of the completely 
longitudinal pole vertex V, which is intimately associated with the Schwinger mechanism of 
gauge-boson mass generation. In that sense, the situation is completely analogous to what 
happens with standard observables, where all direct effects from this particular vertex vanish, 
due to current conservation or general on-shellness conditions. Note that the same situation 
applies to the case of conventional three-gluon vertex, and its associated V; they too cancel 
out completely from the corresponding lattice quantity R. Therefore, any (apparently) 
singular behavior that may be observed in simulations of these quantities should not be 
interpreted as a potential consequence of the pole vertices. 

A possible determination of J m (g 2 ) can provide considerable theoretical insights that 
extend beyond the accurate extraction of the effective charge, and could help us explore 
the nonperturbative behavior of additional key dynamical ingredients. In particular, the 
complete integral equation that govern the evolution of J m (q 2 ) depends on the fully dressed 
four-gluon vertex that involves one background and three quantum gluons. The available 
information on this vertex is very limited at the moment; the only robust result known is 
the WI that is satisfies when contracted with the momentum of the background gluon. It is 
therefore reasonable to expect that further research will be devoted in this direction. Then 
it is clear that independent information on J m {<l 2 ) m ay prove valuable for determining the 
structure of this elusive vertex, at least in some simple kinematic limits. 

In general, the lattice simulation of the PT-BFM propagator and vertices would offer the 
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unique opportunity to verify explicitly powerful formal relations, connecting basic Green's 
function of the theory. Note, for instance, that the simulation of A(g 2 ) could furnish a 
direct confirmation of the fundamental BQI of Eq. (12.91) . given that both A(g 2 ) and G(q 2 ) 
have already been simulated on the lattice. In addition, the result of Eq. (I5.37P may be 
of certain usefulness for future lattice endeavors. Specifically, one may use combinations 
of lattice results to probe the veracity of the (few) theoretical assumptions entering in its 
derivation; conversely, one may assume the validity of Eq. ( 15 . 3 7|) and use it to validate 
the lattice implementation of the BFM algorithm. To be sure, such comparisons must be 
carried out between lattice simulations possessing similar parameters (e.g., bare couplings, 
spacings, volumes). 

Finally, note that the present considerations may be extended to include other fundamen- 
tal vertices of the PT-BFM formalism, such as the vertex connecting a background gluon 
and a ghost-anti- ghost pair, usually denoted by T^. This vertex has a rather reduced tenso- 
rial structure, and satisfies a simple Abelian WI, relating its divergence to the difference of 
two inverse ghost propagators. A preliminary study reveals that relations similar to (I4.15P 
may be also obtained for the ghost dressing function. Hence, lattice simulation of T^ may 
provide nontrivial cross-checks on the infrared behavior of this latter quantity. We hope to 
present the full details of the related analysis in the near future. 
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